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ABSTRACT 

We have studied the variability of 6 low redshift, radio quiet 'PC quasars on 
three timescales (days, weeks, and months) using the ROSAT HRI. The quasars 
were chosen to lie at the two extreme ends of the ROSAT PSPC spectral index 
distribution and hence of the H/3 FWHM distribution. The observation strat- 
egy has been carefully designed to provide even sampling on these three basic 
timescales and to provide a uniform sampling among the quasars 

We have found clear evidence that the X-ray steep, narrow H/3, quasars 
systematically show larger amplitude variations than the X-ray flat broad H/3 
quasars on timescales from 2 days to 20 days. On longer timescales we do not 
find significant differences between steep and flat quasars, although the statistics 
are poorer. We suggest that the above correlation between variability properties 
and spectral steepness can be explained in a scenario in which the X-ray steep, 
narrow line objects are in a higher L/L Edd state with respect to the X-ray flat, 
broad line objects. 

We evaluated the power spectrum of PG1440+356 (the brigthest quasar in 
our sample) between 2 x 10~ 7 and ~ 10~ 3 Hz, where it goes into the noise. The 
power spectrum is roughly consistent with a 1/f law between 10~ 3 and 2 x 10~ 6 
Hz. Below this frequency it flattens significantly. 



Subject headings: quasars — variability, X-rays, emission lines 



- 2 - 



1. Introduction 

While X-ray variability in Seyfert galaxies has been the subject of intensive study (see 
review by Mushotzky, Done & Pounds 1993, and Green et al 1993, Nandra et al. 1997), 
the variability of quasars, being fainter, have received far less attention (the Zamorani et al. 
1984 study remains the most extensive to date). Higher luminosity AGN had been expected 
to be physically larger and so have slower, and most likely lower amplitude, variations. 
A typical quasar is 100-1000 times more luminous than the highly variable Seyferts, such 
as NGC 4051, and so would vary at a significantly lower rate, i.e. a minimum of several 
days. It was probably this expectation that deterred extensive observing campaigns. This 
is unfortunate since we show here that X-ray variability is common, rapid and of quite large 
amplitude (> factor of 2). The variability most likely originates in the innermost regions 
of quasar, and so can help unravel the basic parameters of the quasar central engine (mass, 
geometry, radiation mechanisms, radiative transfer) none of which are yet well constrained. 

These speculations are supported by compact galactic sources, for which the investi- 
gation of the "variability properties vs. spectral shape" and the "variability properties vs. 
luminosity" planes have brought a great improvement in our understanding (see e.g. van 
der Klis 1995). The same is likely to happen for quasars. Compact galactic sources are 
usually bright and variable on timescales as short as 1 ms, which means that their variability 
timescales can be probed by just a few observations of individual objects. For example, the 
Galactic black hole candidate (BHC) Cyg X-l recently underwent two dramatic changes in 
its variability and spectral properties (see e.g. Cui et al. 1997), and it was possible to follow 
the whole cycle from the usual low and hard state, to a medium-high, softer state, and back 
to the low and hard state in the course of a few months. This is unlikely to happen in 
a luminous AGN, even if there were similarities between AGN and black hole candidates, 
because the variability timescales (and the sizes and luminosities) are probably much longer 
(greater and higher) in AGN. Instead, if the analogy between quasars and BHC holds, we 
will observe a population of quasars some in a 'high and soft' state, and some in a 'low and 
hard' state. 

Hence the analogous observational way forward in quasars is to investigate the vari- 
ability properties of samples of quasars, selected according to their spectral properties and 
luminosity. 

The evidence for rapid, large amplitude, X-ray variability in a few AGN with unusu- 
ally steep X-ray spectra and unusually narrow Balmer lines (mostly Narrow Line Seyfert 1 
galaxies, NLSyl, given the strong correlation between these two quantities found by Laor et 
al., 1994, 1997, Boiler, Brandt & Fink 1996) has recently grown: 
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1. NGC4051 shows large variations (50 %) on timescales of ss 100 seconds and has very- 
narrow optical and UV emission lines and steep 0.1-2 keV X-ray spectrum. It is also 
highly variable in the EUV (a factor of factor of 20 in 8 hours, Fruscione et al 1998). 

2. IRAS13224-3809 has a particularly steep 0.1-2 keV spectrum (a x > 3) narrow H/3 
line, very strong Fell emission and shows X-ray variations of a factor 50 or more on 
timescales of a few days (Otani 1995, Brandt et al. 1995, Boiler et al 1997) and a 
factor of 2 in less than 800 seconds (Boiler et al. 1996); 

3. RE J 1237+264 showed one very large (factor of 50) variability event (Brandt, Pounds 
& Fink, 1995) 

4. PHL1092 varied by a factor of 4 in 2 days (Forster & Halpern 1996, Lawrence et al. 
1997). PHL1092 is a high luminosity (5 x 10 46 erg s -1 ) very steep 0.1-2 keV spectrum 
and narrow emission line quasar. 

5. The relatively high luminosity quasar NAB0205+024 (Lo.5-iofc e y = 8 x 10 44 erg s _1 ) 
shows variations of a factor of 2 in less than 20 ks in an ASCA observation (Fiore et 
al. 1998a). 

6. Mark 478 (PG1440+356) shows factor of 7 variations in 1 day during a long EUVE 
monitoring (Marshall et al 1996) 

7. The extremely soft NLSyl WPVS007 (ao.i-2fc e v = 7.3) showed a huge variation (fac- 
tor of 400) variation between the RASS observation and a follow-up PSPC pointed 
observation taken 2 years later (Groupe et al. 1995). 

At this point we believe that a systematic study of a sample of normal quasars spaning 
the observed range of properties is needed. A systematic study of the variability properties 
of AGNs is not an easy task, since it requires: 

1. the selection of well defined and representative (i.e. unbiased) sample; 

2. the availability of numerous repeated observations; 

3. a carefully designed observational strategy. In particular, it is very important that 
the sampling time is regular and that it is similar for all objects in the sample, so 
that the results on each object can be easily compared with each other. Otherwise the 
differences in the observed variability may be induced just by differences in sampling 
patterns. 
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Because of these stringent requirements systematic studies are still lacking (for example 
the Zamorani et al. 1984 study on quasars, and the Green et al. 1993 study on Seyfert 
galaxies do not fulfill any of the three criteria above). To explore in a systematic manner the 
possible relation between line width, soft X-ray spectrum and X-ray variability properties, 
we initiated a pilot program using the ROSAT HRI. In the following we present the results 
from a campaign of observations of six PG quasars which addresses the above points. 

2. The sample 

Six quasars were chosen from a complete sample of 23 optically selected (PG) quasars 
(M B < -23), all studied with the ROSAT PSPC (Laor et al. 1994, 1997). To maximize 
the rest frame soft X-ray detectability in the whole 0.1-2 keV PSPC band the Laor et al. 
sample of PG quasars (M B < —23) is limited to quasars with: low redshift (z<0.4) and low 
Galactic N# (< 1.9 x 10 20 cm -2 ). This sample of AGN has the advantage that it is extracted 
from a complete sample of optically selected quasars, and should therefore be representative 
of such quasars. In particular, the Laor et al. sample selection criteria are independent of 
the X-ray properties, and thus the sample is free from X-ray selection biases. 

The six quasars were chosen to lie at the two extreme ends of the PSPC spectral index 
distribution (/„ oc v~ ax , 1.2<ckx~2.1) and hence of the H/3 FWHM distribution. Figure [3] 
show the H/3 FWHM as a function of ax for the radio quiet quasars in the Laor et al sample. 
(We excluded the X-ray weak quasar PG1001+054 which is too faint for followup variability 
studies with the HRI). X-ray steep and flat quasars discussed in this paper are identified 
by filled squares and open circles respectively. The quasar with ax < 1 and large errors is 
PG1114+445, which has an ionized absorber along the line of sight (Laor et al. 1997, George 
et al. 1997). 

Objects with extreme properties should help highlight the basic correlations which may 
provide important hints for the underlying physical processes. Table 1 gives the redshift, 
the PSPC energy index, the 2 keV and 3000 A luminosities (for Hq = 50 km s~ x Mpc -1 , 
go = 0.5) and the Hp FWHM for the six quasars. All the chosen quasars are radio-quiet. 

3. The observational strategy 

The observational strategy was designed to: a) provide even sampling on three basic 
timescales (days, weeks, and months), b) provide a uniform sampling among the quasars, 
and c) keep the total integration time reasonably small. These constraints result in a quasi- 
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Fig. 1. — H/9 FWHM as a function of ax for the quasars in the Laor et al sample. Filled 
squares and open circles identify the three X-ray steep and the three X-ray flat quasars 
discussed in this paper. 
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Table 1: Continuum and Line Parameters 



name 


z 


a x 


vL v (2keV) 


i/L„(3000A) 


H/5FWHM 








erg s _1 




erg s 1 


km s _1 






Steep X-ray quasars 






PG1115+407 


0.154 


1.9 


43.61 




44.96 


1720 


PG1402+261 


0.164 


1.9 


44.08 




45.45 


1910 


PG1440+356 


0.077 


2.1 


43.64 




44.91 


1450 






Flat X-ray quasars 






PG1048+342 


0.167 


1.4 


43.85 




45.14 


3600 


PG1202+281 


0.165 


1.2 


44.25 




45.00 


5050 


PG1216+069 


0.334 


1.4 


44.78 




45.94 


5190 



logarithmic observation plan: one ~ 2000 s exposure a day for one week; one ~ 2000 s 
exposure each week for one month; two ~ 2000 s exposures a month apart, 6 months after 
the main campaign, i.e. 6-7 points with an even spacing of 1 day; 4-5 points with an even 
spacing of 1 week 2 points with a spacing of 1 month and 2 sets of 2 points with a spacing 
of 6 months. The S/N obtained allows the 3a detection of ~ 30% variability in the faintest 
source of the sample. 

4. Results 

Tables 2 and 3, respectively, give the observation date, the net HRI exposure and the 
background subtracted HRI count rate for the X-ray steep and flat quasars. Counts were 
extracted from a region of 60 arcsec radius. Since ~ 90% of the counts from a point source 
fall within the same extraction region (David et al. 1997) no correction has been made for 
lost flux. Background was extracted from two 60 arcsec radius regions straddling the source 
region. The background is small in all cases (< 10 counts ks" 1 in a 60 arcsec radius regions). 

We have also looked at longer variability timescales (of the order of 2000 days) by 
comparing the HRI fluxes with the PSPC fluxes obtained during the Rosat All Sky Survey 
(RASS) and the pointed ROSAT PSPC observations reported in Laor et. (1997). We 
converted the PSPC count rates to effective HRI count rates by assuming the spectral shape 
observed in the ROSAT pointed observations (Laor et al. 1997). RASS and pointed PSPC 
observations 'converted' count rates are given in Table 4. 
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4.1. Logarithmic normalized light curves 

Figure shows the logarithm of the light curves of the six quasars, each normalized 
to the mean flux, spanning about one month each, while figure |3| shows the observations 
performed six months later compared with the average flux and the dispersion observed 
on timescales shorther than one month. Figure § shows the normalized RASS and PSPC 
pointed observation points, the dispersion as in Figure |2] and the normalized HRI observation 
performed 6 months later. 

The steep ax quasars show large variations (up to a factor of 2) on all the timescales 
investigated. Conversely, the flat ax quasars show little variability on short timescales. 
Factor of 2 variations with respect to the average flux seem to be allowed for both X-ray 
steep and X-ray flat quasars on 200-2000 days timescales, although the statistics on these 
timescales are poorer. 

4.2. Structure function 

To compare the variability of the quasar groups more quantitatively we computed the 
logarithm of the ratio between each pair of flux measurements (at two times U and tj in the 
quasar rest frame). We then computed the median and the mean in 9 time bins for both flat 
and steep X-ray quasar samples: 

(1) < Am >= median or mean(\2.5 * log(f (tj) / f (t(i)))\) . 

The function using the mean is similar to the so called 'average structure function' (see 
e.g. Di Clemente et al. 1996). 

We have verified that no single quasar dominates the median and mean in eq. (1). 

< Am > for the median is shown in figure |5|. Error bars represent here the semi 
interquartile range. < Am > for the mean is shown in figure |6|. Error bars represent here 
the error on the mean and are therefore less conservative. 

It is clear that the < Am > (both mean and median) for X-ray steep quasars is sig- 
nificantly higher than that of X-ray flat quasars on timescales from 2 days to 20 days. On 
the longest timescale (200 days) the two < Am > are consistent with each other, with large 
errorbars. The point at about 2000 days in figures [| and [| compares the HRI monitoring 
results with the RASS results. The point at about 1000 days compares the HRI monitoring 
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Fig. 2. — The ROSAT HRI light curves (plotted as logarithm of the ratio of the flux to the 
average flux) of the six quasars on a 2-20 day timescale. 
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Fig. 3. — The HRI observations of the six quasars performed six months later compared 
with the average flux and the dispersion observed on timescales shorter than one month 
(20-200 days timescale). 



-10- 



Steep 



j 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

- = PG1440+356 



0.2 



Flat 



i 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

PG1202+281 



A 
V 

\ 
u. 

uo 
o 



-0.2 



-0.4 - 



500 1000 1500 2000 
Day (Jan 1990) 



2500 



-0.4 



500 1000 1500 2000 
Day (Jan 1990) 



2500 



1 I 1 I | I I 1 I | 1 I 1 I | I I I I | I I I I 

PG1 115+407 

F 







500 1000 1500 2000 2500 
Day (Jan 1990) 



0.2 - 



-0.2 



II I | I I [ I | II II | II II | I II 

PG1048+342 







500 1000 1500 2000 2500 
Day (Jan 1990) 



i i i i I i i i i I i i i i I i — i i i | i i i i 
! PG1402+261 ; 



0.2 - 



I I I I | I I I I | I 1 I 1 | I I I I | 1 I I I 

PG1216+069 



v 



-0.2 - 



500 1000 1500 2000 
Day (Jan 1990) 



2500 



500 1000 1500 2000 
Day (Jan 1990) 



2500 



Fig. 4. — The HRI observations of the six quasars compared with the RASS and with the 
PSPC pointed observations (200-2000 days timescale). 
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Fig. 5.— "Median" Structure Function (< Am >, see eq. 1) for the X-ray steep (filled 
squares) and flat (open circles). X-ray flat quasars points have been slightly shifted from 
real At for a sake of clarity. 
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'Mean" Structure Function 
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Fig. 6. — Mean Structure Function (< Am >, see eq. 1) for the X-ray steep (filled squares) 
and flat (open circles). X-ray flat quasars points have been slightly shifted from real At for 
a sake of clarity. 
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and the RASS results with the ROSAT PSPC pointed observations. The error bars here are 
larger because of the different sampling times used. 

4.3. Power spectrum of PG1440+356 

For the brightest quasar in our sample (PG1440+356) it is possible to investigate 
timescales shorter than the duration of each single HRI exposure. Analysis of the HRI 
light curves of each single OBI reveals significant variability down to a time scale of a few 
hundred seconds. To quantify the amplitude of this variability we have computed a power 
spectrum for each of the single HRI exposures and, to improve the statistics, we have added 
together the 12 power spectra. The resulting coadded power spectrum is shown in figure 0, 
along with the lower frequency power spectrum calculated at four reasonably well sampled 
frequencies using the full HRI monitoring. The error bars represent the dispersion of the 
power in each bin. 

The power spectrum above 10 -3 Hz is flat, and thus consistent with white noise. Excess 
power is detected at ~ 5 x 10~ 4 Hz, which is a factor ~ 100 smaller than that detected at 
~ 5 x 10~ 6 Hz, a mean 1/f distribution. Unfortunately we do not have information on the 
power in the decade between 10 -5 and 10~ 4 Hz and therefore we cannot tell whether the 
power spectrum decreases smoothly as 1/f above ~ 5 x 10 -6 Hz or rather if it breaks at an 
intermediate frequency. Continuous monitoring up to several days is needed to answer to 
this question. 

The uncertainty on the power below 10~ 6 Hz is large and then we cannot assess whether 
the power spectrum stays constant below this frequency or rather decreases, as suggested 
by the low frequency point in figure [7|. We can however exclude that the power spectrum 
continues to increase as 1/f below 10~ 6 Hz, although the precise break frequency is ill defined. 
Systematic monitoring on weeks to months timescales can strongly constrain this frequency. 
In the framework of accretion disks models, this timescale may be related with the thermal 
or viscous instabilities, and its measure could help to estimate important model parameters 
like the mass of the central black hole and the viscosity parameter a (e.g. Siemiginowska & 
Czerny 1989) 

5. Comparison with other work 

Nandra et al. (1997) used the source root mean square variation a^ ms from ASCA SIS 
data for a sample of (X-ray bright) AGN to suggest an inverse relation between variability 
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Fig. 7.- The power spectrum of PG1440+356 from 0.03 to 2 x 10" 7 Hz. The power 
spectrum is normalized such that the level of the noise is 2 (dashed line). 
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and X-ray luminosity. Our sample of high luminosity, optically selected objects allows us to 
test this suggestion. We have calculated the noise subtracted of ms according to the method 
of Nandra et al (1997). These are reported in Table 5 for two different time scales: days and 
months-to-years. (We do not show errors on of ms because the analytic formula provided by 
Nandra et al. 1997 is not valid in the limit of small number of points.) The light curve of 
PG1216+069 in ||] is consistent with a constant value and therefore of ms is formally zero in 
this case. We take the value of of ms calculated for the other flat X-ray quasars (~ 0.01) as 
an upper limit of of ms in PG1216+069. 

Figure [8] shows of ms as a function of the X-ray luminosity and of the soft (0.1-2 keV) 
X-ray spectral index for the six quasars in our study and for the sources in Nandra et al 
(1997) with no strong low energy absorption (to provide a low energy X-ray spectral index). 
For this figure we used the of ms SIS 0.5-2 keV in Nandra et al (1997). This band overlaps 
quite well with the HRI band (the HRI has some response in the 'carbon' band, below 0.3 
keV, but this is much smaller than in the 0.5-2 keV band). The PSPC spectral indices 
for the Nandra et al. objects are taken from Walter & Fink (1993) and Laor et al (1997) 
or are calculated from WGACAT hardness ratios when not available in these papers. The 
variability of the Nandra et al. objects was measured on a time scale of ~ 1 — 4 x 10 4 s, i.e. 
about an order of magnitude shorter than the timescales probed for the PG quasars (for the 
six PG quasars a^ ms in figured refers to the 2-20 days time scale). 

The addition of the six quasars makes a substantial difference. For the Nandra et al. 
sample alone the of ms vs L\k e v linear correlation coefficient is r=-0.82 (probability of 99.9 
%), while with the quasars it is reduced to r=-0.54. This is still significant (probability 
of 95 %), although it depends entirely on one point (the low luminosity Seyfert 1 galaxy 
NGC4051). The values of cil ms and ax were uncorrelated for the Nandra et al. sample alone 
(r=+0.37, probability of 15%); adding the quasars results in r=+0.55, a 95 % probability 
correlation. 

At a given luminosity the three x-ray 'steep' quasars from our sample all have of ms 
higher than both the 'flat' X-ray quasars in our sample and all the other objects in the 
Nandra et al. (1997) sample by a factor of 10 to 100. The u 2 rms for the Nandra et al. (1997) 
sample refers to a factor of ten smaller timescale, and this might explain why they show lower 
amplitude variability for the same X-ray luminosity. However, if the power spectrum of the 
steep X-ray spectrum quasars falls like « 1// on the 0.1 to 1 day timescales (as suggested 
for PG1440+356, see §4.3), then only part of the difference in a^ ms could be explained. If 
the power spectrum stays more or less constant on these timescales then this difference is 
highly significant. A comparison of variability amplitude in this frequency range between 
broad and narrow line Seyfert galaxies using ASCA and BeppoSAX light curves of similar 
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Fig. 8. — The 'excess' variance as a function of the 1 keV luminosity (in erg s _1 Hz -1 ) 
and of the 0.1-2 PSPC spectral index. Open triangles show quasars and Seyfert galaxies 
from Nandra et al. 1997, filled squares identify steep X-ray spectrum quasars, open circles 
identify flat X-ray spectrum quasars. 
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sampling and length is in progress and will be part of a forthcoming paper (Fiore et al. 1998, 
in preparation). 

Finally we note that the 1 keV X-ray luminosity and the spectral index are slightly 
anti-correlated (r=-0.20) in the sense that steeper objects tend to have lower 1 keV lumi- 
nosities. The correlation is not significant here, but Laor et al. (1997) find a more significant 
correlation between ax and the 2 keV luminosity in their complete, optically selected sample 
of PG quasars (r s = 0.482, probability of 98 %). Thus, low luminosity AGNs tend to be 
more variable because of the luminosity-variability correlation, but on the other hand they 
also tend to have flatter spectra, which diminishes their variability due to the spectral index- 
variability correlation. To remove the effect of the ax — L x correlation from the of ms — L x 
and ol m -ax correlations we performed a partial correlation analysis (e.g., Kendall & Stuart 
1979). The correlation coefficients turns out to be very similar to the previous ones (-0.53 
and 0.54), giving a probability for a correlation again of 95 %. This is still a marginal result, 
which could be highly strengthened or weakened by adding just a few more points at low 
and high luminosities and at flat and steep «j. 

A correlation between variability and spectral shape, similar to that discussed in this 
paper, has been reported also by Green et al. (1993) and Konig, Staubert and Timmer 
(1997) in two different analyses of EXOSAT lightcurves of AGN. 

Di Clemente et al. (1996) studied the variability of 30 PG quasars in the red band and 
compared it with that of a sample of 21 PG quasars observed by IUE. Their full sample 
contains 39 quasars, the majority of which (24) at z < 0.4, as the quasars discussed in this 
paper. Three of the six quasars discussed in this paper are also part of the Di Clemente et al 
sample (the three X-ray flat quasars). Di Clemente et al. found that the variability amplitude 
increases with the rest frame frequency on both 0.3 and 2 years timescales. Extrapolating 
their correlation in the X-ray band one would predict a < m > of 0.75 and 0.83 for the 2 
timescales. This is somewhat surprisingly close to the measured < m > on these timescales, 
given that variability in the red band, UV and X-ray may well have a completely different 
origin. 

6. Discussion 

In our sample of six PG quasars we found clear evidence that X-ray steep quasars 
show larger amplitude variations than X-ray flat quasars on timescales from 2 days to 20 
days. On longer timescales we do not find significant differences between steep and flat 
quasars, although the statistics are poorer. While the sample in this pilot study is small 
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(and expanded monitoring of a larger sample of quasars is surely needed), the distinction is 
so clean cut that we feel justified in speculating on its origin. 

Large narrow band flux variability in sources with a steep spectrum could be induced by 
small changes in the spectral index without large variations of the spectrum normalization. 
If this is the case each source should show large spectral variability with a x anti-correlated 
(correlated) with the observed flux if the pivot is at energies lower (higher) than the observed 
band. Spectral variability of this kind has not been observed in the PSPC observations of 
these and other quasars (Laor et al. 1997, Fiore et al. 1994) and in NLSyl (Boiler et al 1997, 
Fiore et al. 1998a, Brandt et al. 1995). Although our HRI observations cannot directly rule 
out this possibility, we therefore conclude that the observed temporal variability pattern is 
not likely the result of complex spectral variability. 

Laor et al (1997) suggest that a possible explanation for the remarkably strong ax-Hp 
FWHM correlation is a dependence of ax on L/L^dd- The line width is inversely proportional 
to y/I/L Edd if the broad line region is virialized and if its size is determined by the central 
source luminosity (see Laor et al. 1997, §4.7). So narrow-line, steep (0.1-2 keV) spectrum 
AGNs emit close to the Eddington luminosity and have a relatively low mass black hole. A 
similar dependence of spectral shape on L/Lvdd is seen in Galactic black hole candidates 
(BHC) as they change from the 'soft-high' state to the 'hard-low' state. 

A physical interpretation for this effect, as described by Pounds et al. (1995), is that 
the hard X-ray power-law is produced by Comptonization in a hot corona and that as the 
object becomes more luminous in the optical-UV, Compton cooling of the corona increases, 
the corona becomes colder, thus producing a steeper X-ray power-law. In BHC in 'soft-high' 
states this power law component emerges above ~ 10 keV, while the spectrum below this 
energy is dominated by softer emission often associated with optically thick emission from 
an accretion disk. In this model quasar disc emission is at a too low an energy to observe, 
since the disc temperature scales with the mass of the compact object as M B ]( A , leaving 
the Comptonized power law to dominate the 2-10 keV spectrum. The 0.2-2 keV quasar 
emission could be due to Comptonization (see e.g. Czerny & Elvis 1987 and Fiore et al. 
1995) by a second cooler gas component. The 0.2-2 keV and 2-10 keV spectral indices might 
well be correlated with each other, if the emitting regions are connected or if the emission 
mechanisms know about each other. We have undertaken a campaign of observations of the 
quasars in the Laor et al. (1997) sample with ASCA and BeppoSAX to clarify this point. 
Preliminary results (Fiore et al. 1998b) shows that this is indeed the case: steep ax(PSPC) 
quasars tends to have a steeper hard (2-10 keV) X-ray power-law (although the spread of 
the 2-10 keV indices seems smaller than that of the PSPC indices (as also found by Brandt 
et al. 1997). 
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In a sample of quasars with similar luminosities, those emitting closer to the Eddington 
luminosity will also be those with the smaller black hole and hence smaller X-ray emission 
region. Thus light travel time effects would smear intrinsic X-ray variability up to shorter 
time scales in high L/LEdd objects, compared with low L/LEdd objects. Based on this 
interpretation, and on figures [5] and |6| it appears that the emission region of the steep soft 
X-ray quasars is a factor of ~ 10 smaller than that of flat soft X-ray quasars (about « 10 16 
cm and 10 17 cm respectively). Alternatively, the higher variability of steep spectrum objects 
may be a true intrinsic property, which could be induced by some increased instability in high 
L/LEdd objects. The range of luminosity in our sample is too small to tell if the variability 
amplitude appears to be better correlated with black hole mass, as expected for light travel 
time smearing, or with L/LEdd, which would indicate an intrinsic mechanism. 

It is interesting to note that a completely different analysis, based on the interpretation 
of the optical to X-ray spectral energy distribution in terms of emission from accretion disks 
also suggest a small black hole mass and an high arretion rate in two NLSyl (Siemiginowska 
et al. 1998). 

It is also interesting to note that Cyg X-l in the "high and soft" state (Cui et al. 1997) 
shows a total root mean square variability higher than that measured during periods of 
transitions and in the "low and hard" state. This is due to strong 1/f noise, extending down 
to at least a few 10~ 3 Hz, when the source is the "high and soft" state (Cui et al. 1997). 
When the source is in the "low and hard" state this 1/f noise is not present (see e.g. the 
review of van der Klis 1995). 

Ebisawa (1991) found in a systematic study of Ginga observations of 6 BHC that the time 
scales of variability for the soft and hard components are often different. The soft component 
is usually roughly stable on time scales of 1 day or less, while the hard component exhibits 
large variations down to msec time scales. If the time scales of BHC and quasars scale with 
the mass of the compact object the above two time scales translate to 10 4 years and 0.1 
day respectively for our quasar sample. For a sample of quasars with similar redshifts and 
luminosities this predicts a rather small scatter in the soft component and a bigger scatter in 
the hard component. ROSAT results go in this direction. Laor et al. (1994, 1997) find that 
(for their sample of 23 low-z PG quasars) the scatter in the normalized 2 keV luminosity is 
significantly larger than that in the 0.3 keV luminosity. 

We conclude that the analogy between AGN and Galactic BHC seems to hold qualita- 
tively for their X-ray variability properties. 

An alternative and intriguing possibility to explain the correlation between X-ray vari- 
ability amplitude and spectral shape is that a component generated closer to the black hole 
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dominates the emission of steep ax quasars, as in the spherically converging optically thick 
flow proposed by Chakrabarti and Titarchuk (1995) to characterize BHC in the high and 
soft state. If this is the case then we would again expect that the spectrum of the steep ax 
quasars remain steep above 2 keV (and up to m e c 2 according to Chakrabarti and Titarchuk). 
Observations with the ASCA and BeppoSAX satellites instruments, which are sensitive up 
to 10 keV, (Brandt et al 1997, Fiore et al 1998b, Comastri et al 1998) suggest that this is 
indeed the case. 

Schwartz and Tucker (1988) suggested that AGN emission can be produced by an en- 
semble of acceleration sites (i.e. shock waves) with different electron spectral indices and 
therefore emitting power laws with different photon indices. In this picture the spectrum is 
a quadratic function in log E and the mean index at a given energy arises from the greatest 
number of individual acceleration sites. Variability would be greatest for both steepest en- 
ergy index and flattest energy index sources, each of which are dominated by fewer individual 
regions. This is contradicted by the present observations, unless the Soft X-ray flat quasars 
become still flatter (and more variable) above 2 keV. High energy X-ray spectroscopy and 
variability studies are again needed to obtain a definitive answer. A RossiXTE monitoring 
campaign of 4 PG quasars, with a sampling similar to that used in the HRI campaign, is in 
progress and could allow us to clarify this point. 



7. Summary 

We have studied the variability of 6 low redshift, radio quiet 'PG' quasars on timescales 
from a few thousand seconds to years using the ROSAT HRI and archive ROSAT PSPC 
observations. The quasars were chosen to lie at the two extreme ends of the ROSAT PSPC 
spectral index distribution and hence of the B.f3 FWHM distribution. The observation strat- 
egy of the ROSAT HRI compaign has been carefully designed to provide even sampling on 
these three basic timescales and to provide a uniform sampling among the quasars 

We have found clear evidence that the X-ray steep, narrow B./3, quasars show larger 
amplitude variations than the X-ray flat broad H/3 quasars on timescales from 2 days to 
20 days. On longer timescales we do not find significant differences between steep and flat 
quasars, although the statistics are poorer. We suggest that the above correlation between 
variability properties and spectral steepness can be explained in a scenario in which the X- 
ray steep, narrow optical line objects are in a higher L/LEdd state with respect to the X-ray 
flat, broad optical line objects. 

Variability on short timescale (a few hundred seconds) is also evident in the steep ax 
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quasars. The HRI sensitivity allows to study the variability of the brightest of these objects 
down to a few hundred seconds. To investigate shorter timescale variability we will have to 
await missions with high throughput like AXAF and XMM. 

The magnitude of the variability of steep a x quasar does not increase on timescales 
longer than 10 days. In particular for PG1440+356 the power spectrum roughly follows a 
1/f law between 10~ 3 and 10~ 6 Hz, and below this frequency flattens significantly indicating 
a break in the power spectrum. The frequency of the break is however ill defined, because 
of the insufficient sampling on timescales between 1 week and several months. 

The fact that flat a x quasars have comparable variability to steep a x quasars at / ~ 1CT 7 
Hz, but significantly smaller variability on / > 10~ 5 Hz suggests that the variability power 
spectrum of flat a x quasars may be similar in shape to that of steep a x quasars, but with a 
break point scaled to lower frequencies. The other option, that the power spectrum in flat 
a x quasars is just scaled down in amplitude for all /, is less favoured by our results. 

A systematic study of a significantly larger sample of quasars is called for. Such a study 
would allow one to test the strength and significance of the a — a x correlation, to establish 
if the differences in variability power spectra is due to a scaling of the break frequency, or of 
the overall variability amplitude, and to eventually understand if the a — alpha x correlation 
is driven primarily by L/L Edd or by the black hole mass. The ideal instrument for such a 
study would be the ROSAT HRI, which is the only X-ray telescope which allows the required 
large number of short exposures. 
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edge support from ADP grant NAG5-3066. A.L. acknowledges support by the fund for the 
promotion of research at the Technion. We thank Lev Titarchuk and Dan Schwartz for useful 
discussions and comments. 
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Table 2: HRI Monitoring of 3 X-ray Steep Quasars 
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20-21/ Jun/96 


3.86 


201 ±8 




03/Jan/96 


3.69 


1339 ± 19 




04/Jan/96 


0.56 


770 ± 40 




05/Jan/96 


2.66 


847 ± 18 




06/Jan/96 


4.02 


880 ± 15 


PG 1440+356 


07/Jan/96 


4.35 


900 ± 15 




08/Jan/96 


3.99 


596 ± 13 




10/Jan/96 


5.51 


1304 ± 16 




16/Jan/96 


3.51 


800 ± 15 




23/Jan/96 


3.50 


561 ± 13 




20/Jun/96 


2.72 


1276 ± 22 


a in ks; b in ks 1 
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Table 3: HRI Monitoring of 3 X-ray Flat Quasars 



Source 


Date Obs 


Exposure" 


HRI Count-rate b 




07/Nov/95 


4.13 


105 + 5 




08/Nov/95 


2.08 


115 + 8 




09/Nov/95 


2.06 


125 + 8 




10/Nov/95 


2.14 


106 + 8 


PG 1048+342 


ll/Nov/95 


3.30 


103 + 6 




14/Nov/95 


3.40 


78 + 5 




21/Nov/95 


1.86 


91 + 8 




28-29/Nov/95 


3.02 


103 + 6 




04/May/96 


4.31 


87 + 5 




02/Jun/96 


2.54 


43 + 5 




25/Nov/95 


2.59 


153 + 8 




02/Dec/95 


3.51 


173 + 8 




03/Dec/95 


2.19 


167+10 




03-04/Dec/95 


2.37 


164 + 9 


PG 1202+281 


04-05/Dec/95 


2.32 


151 + 9 




06/Dec/95 


2.09 


149 + 9 




07/Dec/95 


2.33 


144 + 9 




08/Dec/95 


2.84 


156 + 8 




09/Dec/95 


3.03 


140 + 7 




16/Dcc/95 


2.10 


197+11 




23/May/96 


3.50 


99 + 6 




21/Jun/96 


4.56 


62 + 4 




09/Dec/95 


2.46 


67 + 6 




10/Dec/95 


2.33 


74 + 6 


PG 1216+069 


ll/Dec/95 


2.65 


76 + 6 




12/Dec/95 


3.16 


67 + 5 




17/Dec/95 


3.80 


69 + 5 




23/Dec/95 


0.88 


67+10 




06/Jul/96 


3.79 


80 + 5 



a in ks; b in ks 1 
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Table 4: RASS and PSPC pointed observation 



Source 


RASS 


PSPC pointed 


Count-rate a 


Date Obs 


Count-rate a 


PG1115+407 


108±9 


65±5 


PG1402+261 


133±10 


195±4 


PG1440+356 


1028±18 


367±6 


PG1048+342 


56±6 


51±3 


PG1202+281 


152±10 


107±3 


PG1216+069 


96±8 


99±4 



a 'converted' (see text) PSPC count rates in ks 1 



Table 5: Source root mean square variation 
name ^(2-20 days) a 2 rms (200-2000 days) 



10^ IP" 2 

Steep X-ray quasars 

PG1115+407 3.2 13.8 

PG1402+261 7.3 19.1 

PG1440+356 8.3 15.2 

Flat X-Ray quasars 

PG1048+342 1.2 13.9 

PG1202+281 0.7 9.7 

PG1216+069 0.0 1.8 



